Gate-induced band ferromagnetism in an organic polymer 
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We propose that a chain of five-membered rings (polyaminotriazole) should be ferromagnetic 
with an appropriate doping that is envisaged to be feasible with an FET structure. The ferromag- 
netism is confirmed by a spin density functional calculation, which also shows that ferromagnetism 
survives the Peierls instability. We explain the magnetism in terms of Mielke and Tasaki's flat-band 
ferromagnetism with the Hubbard model. This opens a new possibility of band ferromagnetism in 
purely organic polymers. 

PACS numbers: 75.10.Lp, 71.20.Rv, 71.10Fd 
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Since the discovery of conducting organic polymers 
jl],D, a great variety of studies have been performed to 
pave a new way to further employ polymers or oligomers 
in realizing various functions, such as field-effect transis- 
tors (FETs) |j) or electroluminescent diodes Q. This is 
most recently highlighted by a series of works by Batlogg 
and coworkers [§],[§, where they have demonstrated that 
some molecular crystals (anthracene etc) and polymers 
(polythiophcnc) can not only be metallized, but even ex- 
hibit superconductivity. 

If a gate-induced superconducting plastic is possible, 
a gate-induced ferromagnetic plastic is of no less inter- 
ests. Ferromagnetism, usually a feature of d or f electron 
systems, taking place in n electron systems has indeed 
been a theoretical and experimental challenge, for which 
an enormous amount of studies have been carried out. 
In particular, purely organic ferromagnets are interest- 
ing. A theoretical proposal was made by Shima and one 
of the present authors, where graphites with super- 
honeycomb structures are considered, but this has to do 
with Lieb's ferrimagnetism g, a kind of antiferromag- 
netism for unbalanced numbers of sublattice sites. Or- 
ganics exhibiting (single-)band ferromagnetism have yet 
to be synthesized. || 

Here we propose a novel possibility for a band ferro- 
magnetism in a purely organic polymer that realizes an- 
other class of flat-band ferromagnetism due to Mielke 
and Tasaki mechanism. |[o| We first find that a chain of 
five-membered rings with a right choice of the functional 
group (polyaminotriazole) has a flat band. When this 
band is made half-filled, where the doping is envisaged 
to be realized in an FET structure, we then show that the 
ground state is ferromagnetic with a spin density func- 
tional calculation, which also shows that the magnetism 
survives the Peierls instability. We finally confirm that 
this really has to do with the Mielke- Tasaki mechanism 
by mapping the 7r-orbital system to a tight-binding model 
with the Hubbard repulsion. This way we end up with a 
gate-induced ferromagnetic plastic. 



The problem of ferromagnetism in the Hubbard model, 
despite a long history since the 1960s 0-H, is still some 
way from a full understanding. One sensible way to study 
the problem is to take some characteristic situation where 
ferromagnetism may be rigorously proved. In fact, Mielke 
and independently Tasaki have proposed some classes of 
models where a flat band exists at the bottom jnj, and 
have proved rigorously for the Hubbard model that the 
ground state is fully spin-polarized when the fiat band 
is half-filled. When the flat band lies in between disper- 
sive bands, the rigorous proof becomes inapplicable, but 
ferromagnetism is still expected when what is called the 
local connectivity condition for the basis vectors in the 
flat band is satisfied. This amounts to a condition that 
adjacent "Wannier" orbitals have to overlap, despite the 
flat dispersion, no matter how they are combined to min- 
imize the orbit size, which is why spins tend to align due 
to Pauli's principle. The case of a middle flat-band is ac- 
tually studied in the context of a model atomic quantum 
wire, fl4|| where the ferromagnetism is shown to be re- 
alized unless the repulsive interaction is below a critical 
value. 

Mielke- Tasaki flat bands are usually constructed from 
interferences between nearest-neighbor and second neigh- 
bor transfers, which requires special lattice structures 
especially in spatial dimensions greater than two. On 
the other hand, in quasi-one-dimensional chains we can 
conceive lattices having flat bands relatively easily fl4) , 
since it is easier to make Wannier orbitals which satisfy 
the connectivity condition in a chain than in a network. 
This is why we first look for some organic polymers. Al- 
though we will have to consider weak three-dimensional 
couplings to estimate the actual T c , [15[ the possibility 
of band ferromagnetism in organic polymers should be 
an intriguing avenue to explore. Here we propose that 
one such system has a flat band. We then show from 
both a band structure calculation and an exact diagonal- 
ization for the Hubbard model that the system should 
be ferromagnetic when the band is made half-filled with 
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a gate-induced doping in an FET (field-effect transistor) 
structure recently developed by Batlogg et al. (||| 

We first note that polymers comprising five-membered 
rings should be promising. This is because we find it 
empirically easier to realize flat bands in chains of odd- 
mcmbcrcd rings. Intuitively, odd-mcmbered rings in- 
corporate frustrations for electron transfers, while even- 
membered (i.e., bipartite) rings have an obvious disad- 
vantage of a tendency toward antiferromagnetism when 
the electron-electron repulsion is turned on, which should 
compete with ferromagnetism. The tight-binding model 
on the connected five-membered rings has indeed disper- 
sionless bands in appropriate, realistic conditions, where 
the eigenstate on the flat band indeed satisfies the local 
connectivity condition (i.e., overlapping orbitals, see Fig. 

0). 

So we start with a search for the case of flat bands by 
scanning various five-membered polymers, i.e., polypyr- 
role, polythiophcnc, polytriazole, polyaminotriazole, etc. 
The band structure is obtained with first principles cal- 
culations within the framework of the generalized gradi- 
ent approximation based on the density functional theory 
(which we call GGA-DFT). [|J| We adopt the exchange- 
correlation functional introduced by Perdew, Bruke and 
Wang |l6| and ultra-soft pseudo-potentials [ pT|JT^ ] in a 
separable form. The wave functions are expanded by 
plane waves up to a cut-off energy of 20.25Ry. As for the 
unit cell, five-membered rings usually alternate their di- 
rections in a chain (see Fig. ||) , so that there are two rings 
in a unit cell with the Brillouin zone folded. We take a 
repeated-chain model along directions perpendicular to 
the chain with a sufficiently large repeat distance. The 
atomic configuration as well as the unit cell size along the 
chain are optimized to minimize the ground state energy 
with the conjugate gradient scheme |l9| . 

It turns out that the flat band is rather hard to realize 
even for five-membered chains, which is not too surpris- 
ing since an odd-mcmbered ring is by no means a suffi- 
cient condition for a flat band. However, attaching some 
functional group to each ring helps since this changes 
site energies and transfer integrals, and we have found 
that among the polymers investigated here polyamino- 
triazole (poly(4-amino l,2,4triazole) to be precise, see 
inset of Fig.||) hits the right condition. Figure || shows 
that the top valence band (with two branches due to the 
above-mentioned band folding) has little dispersion (~ 
O(O.leV)) there. 

To induce the flat-band ferromagnetism, this band has 
to be half-filled, while the band is fully filled when un- 
doped. To assess the possibility of the spin polariza- 
tion in hole-doped polyaminotriazole, we have carried out 
the GGA including the spin degrees of freedom based on 
the spin density functional theory (which we call GGA- 
SDFT) . We have focused on the situation where the 
highest valence band is hole doped to make its upper 
branch empty. Figure S shows the band structure in 



which we have taken a polarized state as the initial state. 
The optimized state remains to be polarized, where the 
splitting between the majority-spin and the minority-spin 
bands is ~ 1 eV. This is the same order of the exchange 
splitting estimated in p0| for the atomic quantum wire. 

@ 

We have to be careful when we dope the band, since 
one-dimensional metals are in general unstable against 
the Peierls distortion, where the electronic energy is low- 
ered at the cost of the lattice distortion. So we have 
checked whether the energy gained due to the spin po- 
larization overcomes the Peierls transition. The left panel 
of Fig. H is the band structure of the Peierls-distortion 
allowed, spin-unpolarized state obtained by the GGA- 
DFT. We can see that the Peierls-splitting at X is negligi- 
bly small. Physically, this should be because cr-electrons 
which form the backbone of the structure have a rigid 
enough bonding that can cope with the Peierls instabil- 
ity here. The total energy of the global minimum in the 
GGA-DFT calculation is higher than that of the polar- 
ized optimized state in the GGA-SDFT calculation by 
-400 meV. 

If we allow the spatial distributions of up and down 
spins to be different in the GGA-SDFT starting from 
an initially unpolarized state, the resulting band struc- 
ture has a wide gap at the Fermi level (right panel of 
Fig.||), which we attribute to an antiferromagnetic gap. 
Although the total energy of the antiferromagnetic state 
is lower than that for the uniformly zero spin density (the 
left panel of Fig. |j), the energy is higher than that of 
the optimized polarized state (Fig.||) by ~50 meV. Thus 
we can conclude that the true ground state in the GGA- 
SDFT is the polarized state, so that we can expect that 
polyaminotriazole becomes, when appropriately doped, 
ferromagnetic at sufficiently low temperature. 

Now we come to the key question of whether the fer- 
romagnetism obtained in the band calculation can be 
identified as the flat-band ferromagnetism for the Hub- 
bard model a Id Mielke-Tasaki. To do so we have first 
to map the 7r-electron system to a tight-binding model 
(see the right panel of Fig. ||), with a Hamiltonian, 

Ho = V,„. I'.'^Cj-t + J2i e i n i in tne standard nota- 
tion. We have estimated the values of transfer integrals 
by first identifying tt bands, then characterizing the wave 
functions as bonding or anti-bonding at Brillouin zone 
center and edges. The site energies are estimated from 
the energy levels of isolated atoms. The obtained fig- 
ures are: <cn — tec — if = 2.5 eV, £nn = 3^4 eV, 
eo — £i = — 1 ~ —2 eV. The band structure for this 
tight-binding model, depicted in the left panel of Fig||, 
indeed reproduces the features of that for 7r-electrons in 
polyaminotriazole (Fig. [j]). If we turn to wave functions 
in Fig. ^ the wave functions at T on the flat band in the 
tight-binding model also reproduce the features of those 
on the flat band in the GGA-DFT result for polyamino- 



triazole. So we can conclude that the tight-binding model 
captures the band structure of polyaminotriazole. 

We can then proceed to the question of whether the 
ground state is spin-polarized in the presence of the Hub- 
bard interaction, Jijj — U^2 i rii-[nii. We can draw a 
phase diagram by noting that a flat band arises when 
a set of equations, 



e — ei 

(s - ei + *nn)/(1 - e) 



(1 - £)0nn - (e - ei) ) - tm, 
-£nn/ (eo — £ ) — ^nn(e — eo)) 



are satisfied, where e is the eigenenergy of the flat band 
and tcN = tec = *o(= 1 here) is assumed for simplicity. 
In Fig. we show the phase diagram against U and e\ 
obtained with an exact diagonalizaion calculation for a 
12-site (2 unit cell) Hubbard model for to = 2.5 eV and 
various values of £nn = 3.0 — 4.0 eV, where eo is chosen to 
satisfy the above equations throughout. We can see that 
we have indeed a ferromagnetic phase unless the repul- 
sion is not too strong (i.e., U <U C with U c = 2 ~ 5 eV). 
The presence of a U c is as expected from the above dis- 
cussion for a flat band lying in between dispersive ones. 
Our preliminary quantum Monte Carlo calculation sug- 
gests that two unit cell is sufficient to roughly determine 
the ferromagnetic phase boundary. 

Let us finally comment on the robustness of the flat- 
band ferromagnetism. First, the one-electron dispersion 
does not have to be exactly flat to realize ferromagnetism, 
as confirmed from a number of studies. |2l],^2| Also, 
electron-electron interactions extending beyond the on- 
site do not necessarily degrade the ferromagnetism. This 
has been shown for the extended Hubbard model where 
the off-site repulsion V does not degrade, or in some sit- 
uation even induces, the ferromagnetism |E3|. As for the 
doping dependence, a numerical calculation [^3 has in- 
dicated that the ferromagnetism survives when the flat 
band is shifted away from the half-filling. Therefore, we 
may expect that the appropriate polymers as exemplified 
by polyaminotriazole should have ferromagnetic instabil- 
ities at low temperatures when the system is sufficiently 
doped in e.g. an FET geometry. In view of a quite recent 
finding M that polythiophene may be doped with this ge- 
ometry to realize superconductivity, we expect the gate- 
induced band ferromagnetism in purely organic polymers 
conceived here should be within experimental feasibility. 
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FIG. 1. The connected five-membered rings, where the 
shading indicates a "Wannier" orbital that satisfies the local 
connectivity condition. 
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FIG. 2. The band structure (left panel) and the optimized 
atomic configuration (right) for the (undoped) polyaminotri- 
azole obtained by the GGA-DFT. The solid (dotted) lines 
represent bands with n (a) character. 
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FIG. 3. The band structure of the doped system with an 
optimized structure obtained with the GGA-SDFT. The solid 
(dotted) lines represent the bands with n (a) character. 
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FIG. 5. The energy band dispersion of the tight-binding 
model with eo = —1.43 eV, £j = —0.5 eV, 
tcN = tec = t{ = 2.5 eV, and £nn = 3.0 eV. To facilitate 
comparison with Fig. ^ we have folded the band to have a 
two-ring unit cell. 




FIG. 6. Bloch wave functions on the flat 7r-electron bands 
in polyaminotriazole obtained with the GGA-DFT (top 
panels) are compared with the corresponding ones in the 
tight-binding model (bottom panels). Left/right panels corre- 
spond to two eigenstates on the flat band at T. White/black 
regions (or circles) represent the sign of the wavefunction, 
while the size of the circles the amplitude. 
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FIG. 4. Band structure for the doped system optimized by 
allowing the Peierls distortion in the GGA-DFT (left panel), 
and that for the antiferromagnetic state obtained with the 
GGA-SDFT (right). 





FIG. 7. Phase diagram for the Hubbard model on a lattice 
depicted in Fig. 5 against U and e\ for various values of £nn 
for to = 2.5 eV. The inset shows the relation between ei and 
eo to realize the flat band. 



